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Abstract 

We provide analytical 0(a s ) results for the three polarized decay structure functions 

H ++ , Hqq and H that describe the decay of a polarized W boson into massive 

quark-antiquark pairs. As an application we consider the decay t — > b + W + in- 
volving the helicity fractions p mm of the W + boson followed by the polarized de- 
cay VF + (t) — > qiq2 described by the polarized decay structure functions H mm . We 
thereby determine the 0(a s ) polar angle decay distribution of the cascade decay pro- 
cess t — > b + W + (— > q\q2)- As a second example we analyze quark mass and off-shell 
effects in the cascade decays H — > W~ + W* + (-^ qiq 2 ) and H — > Z + qq). 
For the decays H — > W~ + W* + (— > cb) and H — > Z + Z*(— > bb) we find substan- 
tial deviations from the mass-zero approximation in particular in the vicinity of the 
threshold region. 



1 Introduction 



The polarization of W ± bosons produced in electroweak production processes is in gen- 
eral highly nontrivial. Therefore, the W bosons produced e.g. in pp(pp) — > + X, 
pp{pp) -> W + W~ + X, e + e~ ->■ W + W~, W + W~ + X and t -»■ b + W + in general have a 
highly nontrivial polarization density matrix. Because of this, there is a rich phenomenol- 
ogy of polarization effects in W production and decay to be explored in present and future 
experiments. For example, one would want to compare the results of polarization mea- 
surements with the predictions of the Standard Model (SM) or models beyond the SM. 

The polarization of the bosons can be probed by decay correlations involving the 
decay products of the polarized W boson. Using such decay correlations, first measure- 
ments of the W polarization in pp — > W + X were reported by the CMS Collabo- 
ration pQ and the ATLAS Collaboration \2\. Measurements of the W ± polarization in 
e + e~ — > W + W~ were published in Refs. pH H]- Finally, results of iy + -polarization mea- 
surements in t — )• b + W + were presented e.g. in Refs. (6J El 13 [10] • Ref. [H] provides a 
survey of SM expectations for the polarization of W bosons in various production channels 
at the LHC. 

In the SM the boson decays into quark or lepton pairs. For unpolarized iy ± -boson 
decays the NLO QCD and electroweak corrections to quark and lepton pair production, 
resp., have been given in Ref. [121 HB]- ^ ne radiative corrections in Ref. [T2l [T3] include 
also quark and lepton-mass effects. To our knowledge the radiative corrections to polarized 
iy ± -boson decays including lepton and quark mass effects have not been done up to now. 

This paper is devoted to the evaluation of the NLO QCD corrections to the decays 
of polarized W ± bosons into massive quark-ant iquark pairs W /± (7) — > q\ q<i where the 
diagonal spin density matrix elements of the W boson can be probed through the polar 
angle decay distribution of the final-state quark pair. We augment our results such that 
they can also be applied to the decay of polarized Z decays into massive quark pairs. In 
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order to provide quick access to the importance of quark mass effects in the decays of the 
W ± and Z bosons we have provided a 0(m^./m^) quark mass expansions of our analytical 
results in a separate paper [H]. In a sequel to the present paper we shall calculate the 
corresponding NLO electroweak corrections to the polarized decay jy + (t) — > ^ + vt [15J. 

In the limit m qi =: m 8 —y our results reduce to rather simple forms which agree with 
previous NLO QCD results extracted from the corresponding calculation of (7*,Z)(f) — y 
qq [16j [T7J HU [191 [20]. Quark mass effects are non- negligible even for on-shell W bosons 
with q 2 = royy for the polarized decay iy + (t) — > cb but become even more important for 
lower values of q 2 as for the decays of off-shell W*^ and Z* bosons as they appear e.g. 
in the recently observed discovery channels H — y W ± W* T and H — y ZZ* of a 126 GeV 
Higgs boson J2TJ [22]. Similarly one needs to retain mass effects in the calculation of 
current-current correlators and their corresponding spectral functions which are needed 
for all values of q 2 . Since there have been claims and counterclaims in the literature as 
to the correctness of known results on radiative corrections to scalar (pseudoscalar) and 
vector (axial-vector) current-current spectral functions, we have compared our unpolarized 
results with previously published spectral function results. 

As an illustration of our general decay analysis we consider the cascade decay process 
t — y b + W + followed by W + —y q\ q 2 where the (helicity frame) diagonal density matrix 
elements of the W + boson resulting from the decay process t — y b + W + have been well 
studied in the literature. We thus provide results on the angular decay distribution for the 
sequential cascade decay t — y b + W + (— y q\ q 2 ) for which we discuss NLO QCD radiative 
corrections in the production process t — y b + iy + (f) and in the decay process iy + (t) — y 
qi q~2- As a second example of much topical interest we take the cascade decay processes 
H —y W~ + W* + (— y qxq-z) and H — y Z + Z*(—y qq) where we discuss quark mass and W* 
and Z* off-shell effects on rates and on angular decay distributions. 

We also briefly comment on the nondiagonal density matrix elements of the W bo- 
son which can be probed by azimuthal correlations in the angular decay distribution. A 
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Figure 1: Feynman diagrams for (a) the Born-term contribution and (b) the one-loop QCD 
contribution to the decay process W + — > qi q 2 

measurement of the azimuthal correlations requires the existence of a preferred transverse 
direction which would be provided e.g. by the transverse polarization direction of the polar- 
ized top quark in the decay t(f) — > b+W + (— > q\ q 2 ). In a similar vein a transverse direction 
can be defined in the large-pr M^-boson production in the process pp(pp) — > W + X. 



2 Born-term results 

Let us consider the quark-ant iquark decay of the SM gauge boson W + 

W + {q)^q l { Pl )q 2 {p 2 ) (1) 

as depicted in Fig. [TJ The LO Born-term amplitude is given by 

M(Born) = M^Born)e,(q) = -i^ij MPi)Y^^v 2 (p 2 ) e,(q), (2) 

where gw is the electroweak coupling constant and the Vij are Kobayashi-Maskawa matrix 
elements (q\ — i; q 2 — j)- We define a reduced matrix element Ai 11 by splitting off the 
common coupling factor — igwVij/y/2 and the factor 1/2 from the chiral projector. The 
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reduced Born-term tensor reads 



H^(Born) = N c £ M^(Born)M ]v (Born) 

quark spins 

= N c Tr ((^ + m 1 ) 7 ^(l - 75 )(^ 2 - m 2 )Y{l - 75)) 

= 8N c {p1p v 2 + P v lP ^-p lP2 g^ + is^ Pla p 2 p). (3) 

The Born-term amplitude ([2]) leads to the LO decay width for an on-shell W + boson with 
q 2 = mly (gly = Ana/ sin 2 ^, Hi = mf/q 2 ), 

nBorn) = \^^ 9 ^\V l3 \ 2 N c H»v(Born){-g, v+ q ^ 



3 8nm 2 w 2 c ^ ^ J \ ^ ' m 2 w J 

= ^g 2 w \V lJ \ 2 N c V\(2-Hi-H2-(Hi-H2) 2 ), (4) 

where \p\ = mwV~X/2, and where A is the value of the Kallen function for the decay 
process, 

A = A(l, Hi, Hi) = 1 + Hi + lA ~ 2/^i - 2// 2 - 2/ii/x 2 - (5) 

The rate expression coincides with the Born-term result in Ref. [12]. 

The subject of this paper are the partial decays from states of the W + boson with 
definite m quantum numbers m = ±1,0, i.e. we are interested in the polarized decay 
structure functions 

H ±± = H tlv e"(±)£* v (±), H 00 = H^e»(0)£* v (0). (6) 

We evaluate the polarized decay functions defined in Eq. ([6]) in the rest frame of the W + 
boson with the z' direction defined by the antiquark g 2 lj The rest frame polarization 
vectors and momenta are thus given by 

e"(±) = -4(0;=Fl,-*,0), <f = ( m W, 0, 0, 



V2 

e"(0) = (0; 0,0,l), P % = (E 2 ;0,0,\ P \), (7) 



■""We have chosen the antiquark direction to define the z' axis in analogy to the antilepton £ + in the 
decay W + — > l + v?. One can equally well choose the quark to define the z' axis. The resulting changes in 
the partial hclicity rate functions will be discussed later on. 



where E 2 = m w (l - p,\ + /i 2 )/2 and \p\ = y/q*J A(l, /ii, /i 2 )/2. 

It proves convenient to bring the rest frame projectors P^_ = e fl (±)e*' / (±) and Pqq = 
£ A '(0)e* I/ (0) into a frame-independent covariant form. One has 

p& = ^(pgr+i-PTiP^), 

IPS? = (8) 

where 

= ^'^W/?, (9) 
iv p 

and where the normalization factor iVp is given by Np = {{p2q) 2 ~ p\q 2 \ m the two-body 
case the normalization factor is reduced to Np = y/q^\p\. The covariant form of the pro- 
jectors are particularly convenient in the NLO tree-graph calculation since the covariantly 
projected integrands in the requisite phase space integrations are Lorentz scalars and can 
thus be handled by the standard covariant methods. 

Using either forms for the projectors (jHJ) or (jSJ), one obtains 

H± ± (Born) = 4N c q 2 (l - Hi - /j 2 ± V\) , H 00 (Born) = 4N c q 2 (l - /ii - /i 2 - A). (10) 

Note that the sum H v = H ++ + H (U : unpolarized transverse) and H L = H 00 (L: 

longitudinal) are fed only by the the parity-even VV and AA current products. The 

difference Hp = H ++ — H (F: forward-backward asymmetric) is fed by the parity-odd 

VA current product. 

At threshold, where q 2 —> (m 1 + m 2 ) 2 , with yZ/Tj" + yfjl^ — >■ 1 and A — > 0, one has 

H (Born) = H 00 (Born) = H ++ (Born) = 8N c mim 2 . All three partial helicity rates 

are equal to one another at threshold. This can be understood from the fact that, at 
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threshold, only the vector current- induced LS amplitude (LS) = (01) survives. This leads 
to the equality of the partial helicity rates using simple Clebsch-Gordan algebra. As we 
shall see in the next section, at threshold one looses the analyzing power of the two-fermion 
decay mode, i.e. the angular decay distribution becomes flat at threshold irrespective of 
the polarization of the W + boson. 

In the massless quark limit \i\ = /i 2 = one has H ++ (Born) = 8N c q 2 ^ and 

Hoo(Born) = H (Born) = as expected from the left-chiral nature of the SM current 

(J2]). The finite mass corrections to the LO helicity structure functions are of O(fii) for H ++ 
and H 0Q , i.e. H ++ (Born) = SN c q 2 (l — jii — [i 2 + . . .) and H 00 (Born) = 4N c q 2 (fii + /i 2 + . . .), 

and of O(fif) for H , i.e. H (Born) = 4iV c g 2 (/zi/i2 + ...). For the sum of the three 

polarized decay functions denoted by H u+L one obtains 

Hu+i(Born) = H (Born) + Hoo(Born) + H ++ (Born) 

= 12N c q 2 (1 - fii - fMi ~ A/3) . (11) 

For the sake of completeness we also define a scalar structure function H tt through 
H tt = H^ u £^(t)e* u (t) where £^(t) = (1; 0,0,0) is the rest-frame time- component (scalar) 
polarization vector of the off-shell W + boson. The corresponding covariant projector on 
the scalar structure function reads 

For the LO scalar structure function one obtains 

H tt (Born) = H s (Born) = AN c q 2 (1 - ^ - /i 2 - A) . (13) 

Note that, at the Born-term level, one has H tt (Born) = H Q0 (Born). H tt (Born) vanishes 
for zero quark masses as expected from current conservation in the mass-zero limit. 
The scalar-longitudinal interference term needed later on is projected by 



such that 

Hot(Born) = H t0 (Born) = -4N c q 2 - /i 2 ) y/\. 



(15) 



3 Angular decay distribution and 

the cascade decay t — > b + W + {— > qi qj 

Consider the rest frame decay of a polarized W + with the normalized diagonal spin density 

matrix elements (p++, poo, P ) where the density matrix elements are given in an unprimed 

coordinate system (x,y,z). Then rotate the spin density matrix around the y axis by an 
angle 9 to a primed coordinate system (x',y,z') such that p' m > m i{0) = Pmm,dl nml (9)dl nm ,(9). 
The angular decay distribution is then determined by the product of the decay probability 
H m i m i for the decay W + (m') — > qi q 2 and the relevant diagonal elements of the spin density 
matrix elements p' m / m i(0), all evaluated in the primed system. 

While a decay analysis in the W + rest system is the optimal choice to probe the density 
matrix elements of the W + boson, the polarization of the W + boson can also be detected 
in other coordinate systems. As an example take the cascade decay t — > b + W + (— > £ + v^). 
When analyzed in the top quark rest system, the polarization of the W + will affect the 

energy spectrum of the final lepton, i.e. leptons from p will be more energetic than those 

from p ++ . 

Returning to the analysis in the W + rest frame we mention that the choice of the z and z' 
axes is a matter of convention and convenience and may be dictated by the physics at hand. 
For example, in the process pp(pp) — > W + + X followed by W + — > l + v several unprimed 
rest frame coordinate systems have been discussed in the literature (Collins-Soper frame, 
recoil frame, target 



the lepton direction 



rame, beam frame) whereas the z' direction is conventionally fixed by 

3 



2 NLO results on W polarization effects in pp — > W + X can be found in Refs. 
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In the example discussed further on (t — > b + W + (— > q\ ^2)) the z direction is fixed by 
the momentum direction of the W + in the top quark rest system (helicity system), and 
the z' direction is determined by the momentum direction of the antiquark qi- 

In addition to the normalized spin density matrix elements p mm with p ++ +p 00 +p = 1, 

it is convenient to work in terms of normalized decay functions H mm = H mm / Y^ m H mm 

such that H++ + Ho® + H = 1. According to what was said before, the normalized decay 

distribution is given by 

W(6) = 2 51 Pmmd mm >{0)d mm ,(9)H m/m/ 

m,m'=0,± 

= 1(1 + cos 2 9) (p ++ + (H ++ + H^) + ^ cos9(p ++ - (H ++ - H^) 

3 » 3 

+ - sin 2 9 (p ++ H 00 + p 00 H ++ + p 00 H— + p— H 00 ) + - cos 2 p 00 H 00 

= I cos 2 [ P++ - 2p 00 + P—){H++ - 2H 00 + H-) 

o 

+ - cos 6{p ++ - {H ++ - H__) 

+ \({p++ + 2p o + P—)(H++ + 2Hoo + H-) - 4p o^oo). (16) 



The distribution ([16]) is a second-degree polynomial in cos 9 and therefore has the form of 
a parabola. Integrating over cos 9 one obtains 

Jw(9)dcos9 = l. (17) 

For unpolarized W + decay one has p = p o = P++ = 1/3 which results in a flat decay 

distribution W(9) = 1/2. Similarly, one obtains a flat decay distribution at threshold 

where H = H o = H ++ = 1/3, i.e. W(9) oc (p + p o + P++)/2 = 1/2 irrespective of 

the polarization of the W boson. 

In the zero quark mass limit and to leading order in a s (where H ++ (Born) = 1 and 
Hoo{Born) = H (Born) = 0) the angular decay distribution f lT6|) reduces to 

W{9) = ^(l + cos^) 2 p ++ + ^(l-cos0) 2 p__ + ^sin 2 ^poo, (18) 



a form quite familiar the analysis of the cascade decay t — » b + W + {— > n + ) [5j El [7J [8j [9] . 

Let us now turn to the a s corrections to the polarized decay functions H mm where we 
include quark mass effects. Surprisingly it turns out that the quark mass corrections to 
the leading NLO term set in linearly and carry rather large coefficients. This has to be 
contrasted with the LO and the NLO unpolarized decay term where the mass corrections 
set in quadratically. In fact, expanding the 0(a s ) polarized decay functions H mm listed in 
Sec. 7 up to 0(y/JIi), one obtains (see also Ref. [H] where the expansion is carried out to 
OOn)) 

H ++ = 8N c q 2 
#00 = 8iV c g 2 



l + l + (7r 2 + 16)y^ + 



+ ^(4-27T 2 V^) + ... 

07T v ' 



H__ = 8N c q 2 



+ ^(l + (7r 2 -16) v ^) + ... 



(19) 



The NLO linear mass corrections are proportional to the antiquark mass m2 and are thus 
maximally asymmetric in the quark massesjf] It is apparent that the NLO linear mass 

terms cancel in the sum H ++ + H 00 + H We mention that the leading order 0(/x°) a s 

contributions can also be extracted from the corresponding calculation of (7*, Z){\) — > qq 
in Refs. [TH [T7J [HI [TH [20] when the quark masses are set to zero in these calculations. 
As concerns the leading order a s contributions, the largest contribution occurs for Hqq 
and amounts to 2a s /(37r) = 2.5% with a s (m 2 y ) = 0.117. The a s corrections can be seen 

to sum up to H ++ + Hqq + H ~ (1 + a s /7r), a result which is well familiar from e + e~ 

annihilation into mass-zero quark pairs. 

The NLO linear mass corrections have rather large coefficients. For example for W + — > 
cb and H ++ , which is the only polarized structure function with a sizeable LO contribution, 



3 When one chooses the z 1 direction along the quark direction (called system I in Ref. P3]), the linear 
mass corrections are proportional to the quark mass m\. As discussed in Ref. |14j the polarized decay 
functions Hl nm in this system are obtained from the present results by the substitution H±±(l,2) — > 
ff£- F (2,l) and H^(l,2) -> Hq (2,1) [TJ] where, using the notation of Ref. [U], the polarized decay 
functions described in this paper are denoted by H^ m (l, 2). 
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the linear mass correction amounts to 155 % (with m& = 4.8 GeV and my/ = 80.399 GeV). 
However, the large mass corrections are tempered when one calculates the normalized 
decay functions H mm which enter the normalized angular decay distribution. In fact, one 
obtains (H ++ + H 00 + = 1) 

H ++ = 1 + ^(-5 + (tt 2 + 16) v ^)+... 

a 



#oo = + ^(4-2^^)+... 



H — = + ^(1 + (tt 2 -16)V^)+... (20) 

where we have used a small a s expansion for the ratios H mm / Hu + l. For W + — > cb the linear 
NLO quark mass effects now amount to only 0(35%) of the leading NLO contribution. 
The reason for the reduction of the linear mass effects is that the largest linear mass 
effect resides in the (unnormalized) polarized decay function H ++ which has a sizeable LO 
contribution. 

The normalized angular decay distribution ( fl6|) can be characterized by the convexity 
parameter (see e.g. Ref. [14] ) 

c f = 7(^% = 7<P++ - 2 Poo + P-)(H ++ - 2H 00 + £__). (21) 
^(cos^)^ 4 

When Cf is negative (positive) the angular decay distribution is described by a downward 
(upward) open parabola. As a second global measure we introduce the forward-backward 
asymmetry of the decay distribution defined by 

_ W{F) - W{B) _ 3 ft ft 

Afb ~ W(F)+W(B) ~ 4 (P++ " p - ){H++ ~ (22) 

where W(F) = W(0 <6<tt/2) and W(B) = W(tt/2 < 6 < vr). If there is an extremum 
of the angular decay distribution in the physical range — 1 < cos# < 1, the extremum is 
given by 

Afb (P++ - p—) (H ++ - #__) 



cos 9 



extr c f (p++ - 2p 00 + P— ) (H ++ - 2H 00 + H— ) 
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(23) 



The three measures are not independent since cos 6* = —ApB/cf. 

extr 

In the small a s expansion and neglecting quark mass effects one has 

Of = ~(l-3p 00 )(l-12^), (24) 

4 07T 

Afb = —(/>++- 6 j^), (25) 

4 OTT 



COS 6 



extr (l-3poo) 6tt ; V ' 



The largest a s correction occurs for the convexity parameter c/. Using a s {my V ) = 0.117 
one finds a 7.5 % reduction of Cf through the radiative corrections, i.e. the radiatively 
corrected angular decay distribution becomes flatter by that amount. This flattening is 
clearly discernible in the plot of the cos 8 distribution of the decay shown in Sec. 8. 

Let us now for illustrative purposes turn to a specific example, namely the cascade decay 
t —7- b + W + (— » qi q 2 )- This process is particularly interesting since the NLO radiative QCD 
corrections factorize into initial- and final-state corrections, i.e. there is no NLO cross talk 
between top quark decay and W decay because of colour conservation [25J. 

The spin density matrix elements of the W + in the decay process t — > b + W + are well 
studied. At LO one has (261 



p ++ (Born) = 0.0007, 



Poo(Born) = ^ = 0.696 -> 0.6887, 

1 + 2x 2 



2x 2 

p(Born) = = 0.304 -> 0.3106, (27) 

1 + 2x 2 

where x = m^jm t . For the numerical values we use the central values of my/ = 80.399 ± 
0.025 GeV and m t = 172.0 ± 0.9 ± 1.3 GeV provided by the Particle Data Group [27]. 
At leading order the density matrix element p ++ is not populated because of angular 



momentum conservation in the two-body decay process. In Eq. (127]1 we have also given 

the NLO QCD results indicated by arrows (cf. Refs. pHJ EHJ EH E] ) fl The correction to 
4 Thc NNLO corrections to the spin density matrix elements of the W + have recently been calculated 
in Ref. 1321. 
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p ++ is very small. The absolute corrections to poo and p amount to 0.73% and 0.66% 

and are thus considerably smaller than the final-state mass-zero corrections to H ++ and 
H 00 given in Eq. (1201). 

If a transverse direction can be specified, one can also probe the non-diagonal spin 
density matrix elements p mm i with m ^ w! . The angular decay distribution is then given 
by [33] 

W(6)= £ An^CW^W^e^' 1 *, (28) 

m,m',m" 

where denotes the azimuthal angle between the production and decay plane. For ml ^ m 
there will be the typical pattern of dispersive and absorptive (or CP violating) contribu- 
tions proportional to cos(m — m') and sin(m — m') <p, respectively. We mention that, if one 
generalizes the above example t — > b + W + {— > q± q 2 ) to the decay of a polarized top quark 
t(t) — > b + W + (— > q± ^2), a production plane can be defined with the help of the transverse 
polarization of the top quark. The corresponding polar and azimuthal distributions are 
given in Refs. [SUED]. 



4 One-loop contributions 

For calculational reasons it is convenient to introduce linear combinations of the diagonal 
helicity structure functions H ++ , H and H 00 given by 

Hi = \{H ++ + #__), H 2 = 1 -{H ++ - H__) t H 3 = 1 -{H ++ + H - 2H 00 ). (29) 

The inverse relations read H±± = Hi ± H 2 and H 00 = Hi — H 3 . Note that the linear 
combinations H 2 and H 3 appear as coefficients of the cos 9 and cos 2 9 contributions in the 
angular decay distribution (fT6j) . 

The one- loop QCD correction to the decay process W + — > qi q 2 is shown in Fig. QJb). 
The vertex correction to the Born-term (V — A) vertex factor 

-^v^-t 1 (») 
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can be written as — i(gw / V^VijAT 1 ^. At NLO one finds 



1 L 



l -^{i - 75 ) + Ar£ = (i + A L )r- 2 



+ Stil- 



ls 



2 ^+15 



(31) 



2 2 ' ^ 2 1 «^ 2 

where, as in the LO case, pi and ]?2 are the four-momentum of the up-type quark and the 
down-type antiquark, respecticely. The UV and IR singular parts reside in the Born-term 
like structure In order to regularize the singularities, we use dimensional regulariza- 
tion with D = A — 2e. The UV singularity is removed by UV renormalization while the 
IR singularity will be cancelled by the corresponding contributions from the tree-graph 
contributions. The form factors are in general complex valued, i.e. they contain absorptive 
parts as can be visualized from Fig. QJb). For the present calculation we only consider 
the diagonal helicity rate functions, and thus we only need the real parts of the one-loop 
contributions. One has 



Re A, 



x 



Re A 



R 



ReB\ 



ReB^ 



ReB 2 L 



ReB 2 R 



-£^r(i + £ )( ^ 

2 2 /il + fl2 ~ (fil-fi2) 
e y/\ 



\2 

In 



1-v 



+ -/=(l-/^i -^2) 



.1 + v 
In (1 - (J~Hi 



+ 3 a/A In 

IT,) 2 



1 + v 



'Hi 



In 



1 + v 



(A*i - H2) In 

+ ReL'\ +4 



a^ Cf ~ 
~a^ Cf 



, \//^2 


VI + d, 


2mi 


"1 - 


2/xi + (//i 


g 2 




VA 


2m 2 
q 2 







In 



1 + v 



(1 - Hi + n 2 ) In 



1 - /ii - /i 2 + (1 + Hi - H2) 



a s 2mi 
An q z 
- Hi - 



x 



a Sri 2m 2 
An q A 



H2 + (1 - Hi + H2) 



In 



(2 + Hi ~ H2) In 



x/A 



In 



(2 - Hi + H2) hi 



'Hi 



+ 1 



'Hi 



+ 1 



iHi, 



+ 1 
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Figure 2: Feynman diagrams for the NLO tree-graph contributions to the decay process 
W + ->■ q x q 2 



x 



1 - 2/i 2 + (/ii - /^ ln / I 



1 + V 



+ (1 + 11!- fx 2 ) In 



//Xi 



+ 1 



(32) 



where CV = (A 7 " 2 — 1)/2N C = 4/3. We have introduced a velocity parameter t> defined by 



(33) 



Per se the velocity parameter has no physical meaning except that it reduces to the usual 
velocity v = yjl — Am 2 jq 2 in the equal mass limit. The function ReL' is given in Ap- 
pendix A. The scale /t in ReA L has been introduced to keep the strong coupling constant 
dimensionless in D = 4 — 2e dimensions. The dependence on ji cancels in the sum of the 
one-loop and tree-graph contributions. The one-loop contributions to the helicity structure 
functions finally read 



H x {loop) = 8N c q 2 (l -m - n 2 )ReA L + 16N c q 2 ^JI^ReA R , 
H 2 (loop) = -8N c q 2 V\ReA L , 
H 3 (loop) = 8N c q 2 \ReA L 

+AN c q 2 \ (mi (Re B\ - Re B 2 L ) + m 2 (Re B l R - Re B 



(34) 
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5 Tree-graph contributions 

In accordance with the Lee-Nauenberg theorem, the IR singularities of the one-loop con- 
tribution will have to cancel against the gluon-emission tree-graph contributions depicted 
in Fig. |2j The result of the full phase-space integration can be expressed in terms of the 
decay rate terms £ , . . . , £±, if (0), S e z (0), l[ (0), Sf (0), and I e (0) listed in Appendix A. Again 
we list our results in terms of the auxiliary expressions Hi, H 2 and H 3 defined in Eq. ff29l) . 
One has 

Hi(tree) = JV[4(1 - »i - ii 2 )D s - 4//i(l + 7^ - n 2 )l[ (0) 

- 2^(1 - 12/ix - 2fi 2 - 5fil + Afnfi 2 + fi 2 2 )S[{0) 

- 2//i(6 + 4/i! - 7fi 2 )£i + 2// 2 (2 + 3^i)£ 2 - 2(1 - 11/zj + /i 2 ) VXl , 



H 2 (tree) 



N 



Ay/XD! + 4(1 - 3//i - 2// 2 - //? + ^1)^(0) 



- 2(2 - /ii + fM 2 - n\ + A*i^ 2 )4 - 8A£ 4 

+ 4V^(1 + 2/xa - + 2^A(2 + ^ + /i 2 )£ 2 



+ (3 + 14^ - 3^i + 3/i 2 ) ((1 - ^JTiY - /i 2 

H 3 (tree) = n[a\D s - 12^(1 + 7/n - A* 2 )/f(0) 

- 6y7/I(l - 12//i - 2/i 2 - 5/x? + 4//i// 2 + nl)S[{Q) 



where 



- 2^(20 + 13//! - 24/i 2 + /x? + /ii/i 2 + 4^|)4 



+ 2/i 2 (4 + 12/ii — fx 2 — 4/ii - /ii/i 2 - /i 2 )4 



— 2(3 — 36/ii — + 8/ii/i 2 — // 2 )vA 



(35) 



iV := a 8 N c C F q 2 /(7rVX), 



(36) 
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D s := (l-fi 1 -fi 2 )(D i + S e z (0))-2VXD 

+- ((i + m - f i 2 )e 1 + (i - n x + fi 2 )£ 2 + Vx) , 

Dj := (l- /Ul - At2 )(^ + /f(0))-2v / AD 

+\ ((1 + A*i - A*2Ki + (1 - Mi + ^2 + VX) . (37) 
We have isolated the IR singular parts in D and D e given by 
D := In ( / ) - 1, 

B ' := ln (7x|7s)" m ' + + ^ Li2(1 - c,+) -I Li2(1 - a - ) (38) 

with a + = (1 — a*i — A*2 + VA)/(1 — A*i — A*2 — \/A) = cC 1 - The IR singularity has been 
regularized by a small but finite gluon mass mc = V-^-Q 2 - Since the one-loop calculation has 
been done using dimensional regularization, one needs to convert the IR divergent piece of 
the tree-graph contribution to the corresponding expression in dimensional regularization 
by using the one-loop relation 

lnA= (^pj (I- 7E + ln(47r)). (39) 

6 Total NLO contribution 

Because of the aforementioned Lee-Nauenberg theorem, the IR singularities cancel when 
adding the one-loop and tree-graph contributions. Using the IR finite quantities 

q 2 

A s ■= D s + ^ReA L 

= ^(1 - fii - /i 2 ) (t A + 2S e z (0)) - V\£ A + X - (\ - pn - a* 2 + ^a) h 

1 3 
+ -(/i! - ii 2 )V\£ b + - ((1 + Ati - Af 2 Ki + (1 - Hi + H2Y2 + >/A) , 

q 2 

Al := Di + 2N BbAl 

= tMt) (tA + 2/f (0)) - sf\£ A + \ (l - in ~ H2 + ^A) £ 3 
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1 3 

+ -AjUx - ^)VX£ B + -[(l + /ix- fj. 2 )h + + /i 2 )t2 + VX) (40) 

(£ai V-b and tj± are listed in Appendix A), the total results read 

H x {a a ) = N[A(l-^-ii 2 )A s - 4 Ml (l + 7/n - /x 2 )/f(0) 

- 2^1(1 - 12/ii - 5/^ - 2/i 2 + 4/!i/i 2 + //i)S'i(O) 



H 2 (a s) 



- 2/i! (6 + 4/i! - 7/i 2 )4 + 2/i 2 (2 + 3/ii)£ 2 



- 8/ii/i 2 4 - 2(1 - ll/ii + /x 2 )VA 



iV 



4^^/ + 4(1 - 3/ii - /i? - 2/i 2 + l4)I £ (Q) 



- 2(2 - /ii - /ii + /12 + Mi/*2)^o - 8A£ 4 

+ 4VX(1 + 2/n - /i 2 )£i + 2^(2 + //i + /i 2 )£ 2 



(41) 



+ (3 + 14 v //i7 - 3/ii + 3/i 2 ) ((1 - y/K) 2 ~ M2 



(42) 



H 3 (a s ) 



n[a\A s - 12/n(l + 7/ii - Ai 2 )/f(0) 

- 6^/^1(1 - 12/ii - 5/i? - 2/i 2 + 4/ii/ia + /i 2 )£f(0) 

- 2/ii(20 + 13/ii + /I? - 24/i 2 + + 4/4)£i 

+ 2/i 2 (4 + 12/ii - 4/ii - /i 2 - /ii/i2 - A* 2 )^2 

+ A (/ii + /i 2 - (/ii - /i 2 ) 2 ) 4 - G"i - /i 2 )Av / A£ B 

- 2(3 - 36/ii - A*i + 8/ii/i 2 - /i 2 ) \/A . 



(43) 

For /ii = /i 2 we agree with our previous NLO QCD results on (7*, Z)(f) — > qq [161 El HI 

USED]. 

A further check can be done by comparing the sum of the partial helicity structure 
functions Hjj + l = H ++ + H Q0 + H = 3Hi — H 3 with the corresponding results [121 [13] . 
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For the unpolarized decay function Hu+l(ch s ) we obtain 



H u+L (a s ) 



N 4(3(1 - //! - // 2 ) - A) A s 



+ 2/ii (2 + fl!+ fll~ I8/X2 + + ^1 

+ 2/i 2 (2 - 18/ii + 4/i^ + /i 2 + /ii/i 2 + A<J) ^2 



((1 - /ii- - A)A - 6/ii/i 2 ) £3 + (/"I - /i2)Av / A£ 



+ 2(1- 5/ii - 5/i 2 - A + 6/^1/12) vA 



(44) 



in full agreement with Ref. |13jrl 

7 High-energy and threshold limit 

Since our results are obtained in analytical form, one can study different limiting cases. In 
the high-energy (or mass-zero) limit one needs to expand the Kallen function up to O(fif). 
One has 

v 7 ! = yfl + fij + ii 2 2 - 2//i - 2/i 2 - 2/i ly u 2 = 1 - \t\ - fj. 2 ~ + C(/i, 3 ). (45) 
The high-energy limit of the decay rate terms are given in Appendix B. One has 



This result has already been used in Sec. 3. 

5 We also find agreement with the final result in Ref. [12] after correcting two typos in Eq. (A. 50) of 
Ref. |12j . namely after removing the denominator factors (l + wi) in two of the Spence functions in (A. 50). 
We thank A. Denner for a communication on these typographical errors. 




(46) 
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At threshold one has ^JJ^l J ^ — > 1 and thus A — > 0. Using the results of Appendix C 
one obtains up to 0(a s ) 

H ++ = H 00 = H = H a -> 8iV c g 2 |v^I^+ 8^^=^^ . (47) 

At threshold, all four 0(a s ) helicity rate functions are equal to one another as is true 
at LO (see the pertinent discussion in Sec. 2). Concerning the on-shell decay of the W + 

involving the polarized decay functions H ++ = H 00 = H one thus has a flat angular 

decay distribution at threshold also at NLO. The Coulomb singularity proportional to 
in Eq. ( 1471) signals that perturbation theory breaks down close to threshold. One 
has to use nonperturbative methods to analyze the region close to threshold similar to the 
analysis of e + e~ —ty,Z—}tt close to threshold discussed in Refs. [3U |35l 136] . 



8 Numerical results for off-shell and on-shell 
polarized decay functions 

In this section we present our numerical NLO results for the three helicity rate functions 
H mm for on-shell and off-shell W bosons. We choose the \[q 1 range to extend from threshold 
y/q 2 = rriiy + m c to the maximal energy \fq% = mt — rrib attainable in the decay t — > b + W + . 
In order to highlight quark mass effects we take the decay channel with the highest quark 
masses, namely the channel W + — > cb proportional to (V^) 2 ~ (0.04) 2 . For the quark 
masses we take the pole masses m t = 172.0 GeV, = 4.8 GeV and m c = 1.5 GeV. We let 
a s run with two-loop accuracy. At q 2 = = 80.385 GeV 2 we have a s = 0.117. 

In Figs. El H] and [5] we display the y/q 2 ' dependence of the Born-term and 0(a s ) helicity 

rate functions H Q0 , H and H ++ for the process W + — > cb. We choose to normalize our 

results to the unpolarized Born-term rate function Hu+L{Bom) given in Eq. ( TlTT) . 

Fig. [3]shows that the ratio Hoo(Born) / Hu + L(Born) rapidly approaches the appropriate 
threshold value of 1/3 at the lower end of the spectrum. The corresponding NLO ratio 
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Figure 3: Energy dependence of the normalized coefficient i?oo/ Hu + L(Born) in the interval 
[mb + m c ,m 4 — m&] at LO (dashed lines) and NLO (solid lines). The dotted vertical line in 
Figs. [3H6] marks the position of an on-shell W boson. 
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Figure 4: Energy dependence of the normalized coefficient H /Hu + i(Born) in the interval 

[nib + rn c , m t — nib] at LO (dashed lines) and NLO (solid lines) 
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Figure 5: Energy dependence of the normalized coefficient H ++ /H u+L (Born) in the interval 
[nib + m c , m t — frib] at LO (dashed lines) and NLO (solid lines) 
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quickly approaches +00 at threshold because of the Coulomb singularity in the a s NLO 
one-loop contribution. Towards the higher end of the \[q* spectrum the two ratios quickly 
reach their respective asymptotic values of zero and 2a s /37r. For the maximal energy 
\[q 1 = m t — mb the results are already close to the high-energy limit. The Born-term 
result approaches zero while the 0(a s ) result stays at a finite value 2a s /3n ~ 0.02 (with 

aslmt — rrib) ~ 0.1). For H {Born) Fig. HJshows that, at maximal energy, the high-energy 

result a s /6ir ~ 0.005 is already obtained with high accuracy while the Born-term result 
again approaches zero. Finally, for the normalized coefficient H ++ (Born) one sees from 
Fig. [5] that the Born-term result approaches the value 1 at maximal energy. 

All three plots show that the approach to the high-energy (or mass- zero) limit is rather 
slow for the a s corrections. In particular one is not close to the asymptotic NLO values 
H—/H u+L (Born) ~ (1 + a s /6ir), H 00 /H u+L (Born) ~ 4a s /67r and H ++ /H u+L (Born) ~ 
a s /67r at the on-shell value \fq I = mw indicated by the dotted vertical lines in Figs. [3j 
H] and O The large NLO mass effects even at the scale yfq 1 = m w have been discussed 
before in Sec. 3 and in Ref. [H] where one can find an 0(//j) expansion of the NLO mass 
effects. 

In Fig. [6] we leave out the Born-term contributions and show the NLO corrections to 
^mm(NLO), divided by the sum of these. It is obvious that, at threshold, the effect of 
the Coulomb singularity drops out in this ratio and all three helicity structure functions 
contribute with a relative factor 1/3. On the other end of the spectrum in Fig. [6] the curves 
start their slow approach to the limiting values 1/6 (for H±±) and 4/6 (for H 00 ). 

In Fig. Owe plot the cos# distribution for W{6). It is quite apparent that the distri- 
bution becomes flatter through the radiative corrections. Numerically one finds 



(-0.809; -0.747) 



(48) 



A FB 



(-0.228; -0.225) 



(49) 



cos 6 



max 



(-0.281; -0.302), 



(50) 
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Figure 6: Energy dependence of the NLO corrections to H (NLO), if o(NLO) and 

# ++ (NLO), divided by the sum H u+L = #__(NLO) + # o(NLO) + # ++ (NLO) in the 
interval [mj, + m c ,m t — m,b} 



24 



20 
10 




Vq 2 = m v 



_J I L 




0(cx s ) result 
Born term result 



-0.5 







0.5 



cos 



Figure 7: Normalized angular decay distribution W(9) = W(9)/W at LO (dashed line) 
and NLO (full line). NLO result contains both initial-state and final-state corrections. 



where the first number gives the Born-term result and the second number gives the 0(a s ) 
result, both with fa ^ 0. The negative value of the convexity parameter Cf means that the 
angular decay distribution is given by a downward-open parabola. 



9 The decays H W~ + W* + {^ qiq 2 ) 
and H -> Z + qq) 

In this section we consider quark mass and off-shell effects in the polar angle distribution 
of the decay I¥* + (t) — > qi<l2 where the off-shell W* + is produced in the Higgs decay 
H — > W~ + W* + . We shall also briefly touch on the subject of the three-body decay H — > 
Z + Z*(— 7- qq). The corresponding leptonic modes have recently been observed at the LHC 
and are therefore adequately dubbed "Higgs discovery channels" [2Tj [22] . Off-shell effects in 
these decays will lead to additional scalar and scalar-longitudinal interference contributions 
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in e.g. the off-shell decay W / * + (t) — > q\q 2 well familiar from neutron beta decay and from 
the semileptonic decay 2° — > S + + ^Tv^ [37J, or from the decay B — > D^*> + t~v t [38] . 
The scalar and scalar-longitudinal interference contributions are quadratic in the quark 
masses and can thus be neglected at the scale q 2 = rriyy. However, for the off-shell decay 
H -)■ W~ + W* + the scale is not set by Tfi^y but by the off-shellness of the H^*^ which 
extends from threshold q 2 = (mi + m-i) 2 (maximal recoil point) to the zero recoil point at 
q 2 = {ttlh — iriw) 2 , i-e. one has (m# = 126 GeV) 

(mi + m 2 ) 2 < q 2 < (m H - m w ) 2 . (51) 

One will therefore have to carefully consider quark mass and W* + off-shell effects in the 
q 2 region close to threshold. 

The differential decay distribution for the decay H — > W~W* + (— > 5152) is given by 

dT _ 9 4 ]v ,2 \Pw\\p\ 1 2 (f)) ( z 9] 

dq 2 dcos9 1024vr3 |Kl21 m 2 HV ^ {q 2 - m 2 w ) 2 + m 2 w T 2 w 3 oS ~ shM{ ^ 1 ] 

(g 2 = SmlyGp/ \/2 = 0.4265) where the polar angle decay distribution reads 

W off - S heii(0) =\{~ 9 m ' + ^f-) ( " 9 VV ' + ^W#mV, (53) 

and where \pw\ = A 1 / 2 (m^-, m^/, q 2 ) j2mu and \p\ = y / ? 7 A 1 / 2 (l, /ii, /i 2 )/2 are the magni- 
tudes of the momentum of the W in the H rest system and the momentum of the quarks in 
the W* + rest system, respectively. In Eq. ( 1521) we have integrated out a trivial azimuthal 
angle dependence. The polarization of the W* + is encoded in the density matrix function 
p^y which in turn is determined from the decay H — > W~W* + . The hadron tensor 
contains the decay dynamics of the decay W* + — > q\q 2 as described in Sec. 3. 

One can separate the spin 1 and spin parts of the propagators in Eq. (1531) by writing^ 



^ + <qi) = ( - r < + rr- - - \)). (54) 



6 In the analysis of Refs. [SHI SO] only the spin 1 piece of the propagator is kept which is adequate for 
the zero lepton mass case. 
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Note that, in the product of the two off-shell propagators in Eq. (1521) . the scalar-longitudi- 
nal interference term acquires an extra minus sign. 

The polar angle decay distribution of a spin 1 boson decaying into a quark pair described 
in Sec. 3 will be augmented by the contribution of a scalar-longitudinal interference term 
and a scalar contribution. One has 

3 

W o ff- s h:el\(0) — ^ X, Pmmd mm ,{6) d mm ,(8) H m > m > 
m,m'=0,± 

-77 (l — ) (PtoHto + PotHot) cos 6 + - (l 2~) PttHtt- (55) 

In the next step we calculate the density matrix elements of the off-shell W* + in the 
decay H — > W'W^l'f) where we sum over the three polarization states of the on-shell 
W~ . In the SM the Higgs particle couples to a pair of W bosons via the metric tensor, i.e. 
the matrix element for H — > W~W + is given by 

M = imwgg»v£w- £ * q u , (56) 

where By/- and e q denote the polarization vectors of the on-shell W~ and the off-shell W* + 
boson, respectively. On squaring and summing over the three spin states of the on-shell 
W~ one obtains the density matrix elements 

2 l_^^^\ £ U m) {m% (57) 



Pmm' = m w -g 1 * + 



The square of the coupling factor g does not appear in Eq. fl57|) since we have taken the 
freedom to absorb g 2 in the overall factor in the rate formula ( |52l) . 

We calculate the density matrix elements p mm i in the Higgs rest frame with the z axis 
along the W* + momentum q = px — py/. Let us collect the relevant expressions for the 
four-momentum and the polarization vectors of the W* + boson. One has 



<f = (q ; 0, 0,\Pw\), qo = ^^(m 2 H + q 2 -m 2 w ), e£(±) = (0; =f1, -i, 
^(0) = -^(\Pw\;0,0,%), e£(t) = X? = 4,(50; 0,0,1^1). (5 



'q K ' \/q v H 
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On evaluating Eq. (157)) one obtains 



P++ 



Pot 



P 



Pto 




\Pw 



(59) 



At threshold (maximal recoil) when q 2 — > (mi + m2) 2 , and for rrii — » 0, the longitudinal and 
scalar contributions poo = Ptt = Pto = ( m /f — m w) 2 l^l 2 become dominant. On the other 

end of the q 2 spectrum (zero recoil), where \pw\ = 0, one finds p ++ = p o = P = m w 

and p tt = pto = 0. 

Since the decay H — > W~W* + is parity-conserving, the transverse density matrix 

elements p ++ and p are identical to each other, i.e. one has p ++ — p = 0. This means 

that there is no parity- violating contribution to the cos 8 coefficient in the (first) spin 1 
part of Eq. (155)) (see Eq. (Tl6l) ). The second cos# contribution in Eq. (155]) does not have 
a parity- violating origin but is a parity-odd effect. It arises from the scalar-longitudinal 
interference contribution with J p properties (0 + , 1 _ ) (VV) and (0", 1 + ) (AA) resulting in 
a parity-odd contribution. 

The polarized decay functions H±± and H Q0 have been calculated before. The LO and 
NLO forms of the additional polarized decay functions H tt and H t o can be found in Sec. 2 
and in Appendix D. For the convenience of the reader we list H a and H t o together with 
their O(pi) mass expansion. One has 



H tt = AN c q 2 (l-^-to- A + H 1 s (a s )) 

= 4N c q 2 f/ii+/i 2 + ... + ^(l8pi + \Sp 2 + l2p l In //i + 12/x 2 In p 2 + ■■■)) (60) 



H t0 = H ot = AN c q 2 (-(p 1 ~p 2 )V\ + H^ t (a s )) 
= -4A c g 2 (p 1 -/i 2 + ... 

+ ^(26/i! - Up 2 - 47r(/i! - p 2 ) + 12/xi ln/Xx + I2p 2 \iap 2 + ...)). (61) 




and 
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On integrating Eq. ( 152|) over cos 9 one obtains the differential q 2 rate which is given by 
dT 9 4 it/ 12 \Pw\\p\ 1 



cig 2 10247T 31 ' 121 m 2 H ^(q 2 ~ m 2 w ) 2 + m 2 w T 2 w 

x \{(p++ + Poo + P-)(H ++ + #oo + ) + 3(1 - \) 2 p u H tt ). (62) 

In the zero quark mass limit rrii — >■ where + i^oo + # = 8N c q 2 and if ti = 0, the 

Born-term rate calculated from Eq. (|62|) can be seen to agree with the result of Refs. [JTJ H2] 
when N c is set to one[] 

In our numerical discussion we again concentrate on the mode H — > W~ + W* + {— > cb) 
in order to highlight quark mass effects even if this mode is suppressed by \V c b\ 2 = (0.04) 2 . 
In Fig. |8] we show the q 2 dependence of the rate. Let us begin our discussion with the 
Born-term contributions. In the threshold region, where the longitudinal W* + dominates, 
the rrii ^ differential rate clearly shows the appropriate threshold behaviour 1\p\j \fq 2 ~ = 
A 1 ' 2 (l, fjLi, 112), i-e. the differential rate vanishes at threshold. This vanishing is not seen for 
the rrii = curve. This can be understood by taking the rrii — V limit of A 1 ^ 2 (l, Hi, ^2) 
keeping q 2 small and fixed with the result A 1 ' 2 (l, ^2) 1- For the q 2 = value of the 
differential rrii = rate one then obtains 



^2 



dT 

dq 2 



= TT^I^ 2 ! 2 ^ 32, T W r2 ^ = 9-554 • 10- 11 GeV- 1 . (63) 
9 2 =0 1024tt 3 3 m 3 H m 2 w (m 2 w + T 2 W ) 



in agreement with Fig. |HJ At higher values of q 2 the difference between the rrii = and 
rrii 7^ Born-term curves becomes smaller and smaller. The radiative corrections are 
largest in the threshold region. Away from the threshold region they amount to over 
10 % and are thus considerably larger than what would result from the simple estimate 
otsl^ ~ 3.7%. We mention that the radiative corrections to the LO rrii — curve in Fig. |H] 
is simply given by multiplying the LO result by (1 + a s /ir). 

Fig. [8] also shows that the 0(a s ) rrii 7^ rate does not go to zero at threshold. This can 
be traced to the presence of the NLO chromodynamic Coulomb singularity at threshold. 



r As pointed out in Ref. [41], the corresponding result in Ref. [43] is too small by a factor of 3/4. 
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Figure 8: Differential rate for the three body decay H — > W~ + W* + {— > cb). The three 
curves correspond to (i) Born term (m; = 0) (dotted line) (ii) Born term (m« 7^ 0) (dashed 
line) and (hi) 0(a s ) with (wij 7^ 0) (full line). 

The Coulomb singularity proportional to A -1 / 2 (see Eq. ( 1C4|) ) is cancelled by the overall 

rate factor \p \ = \fq I \ 1 l 2 /2 resulting in a finite contribution at threshold proportional to 

a s . One can estimate the finite threshold value of the 0(a s ) rate by neglecting terms of 

0(q 2 /m^y) in Eq. ( )62|) whence one can express the finite threshold value in terms of the 

LO rrii = contribution in Eq. ( 1631) . One then obtains 

dT 32tt dT . A ^ 

a<r thresh 3 ag z g 2 =0 

Using a s (g 2 = (4.8 + 1.5) 2 GeV 2 ) = 0.165, one obtains approximate agreement with Fig. [HJ 
As has been emphasized before, perturbation theory cannot be trusted in the threshold 
region and therefore the treatment of the decay W* + — > cb requires a nonperturbative 
treatment including a resummation of the chromodynamic Coulomb singularity. The above 
exercise leading to Eq. ( 164"1) merely serves to check on the consistency of our calculation. 
In Fig. |9]we show a plot of the q 2 dependence of the convexity parameter. The convexity 
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Figure 9: Convexity parameter Cf(q 2 ) as a function of q 2 . Labelling of curves as in Fig. [8] 

parameter is obtained from Eq. (1621) by replacing (p++ + poo + P + -^oo + H ) by 

3/4(p ++ — 2p 00 + p )(H ++ — 2H Q0 + H ), setting the scalar contribution to zero, and 

then dividing by the differential rate (162]) . At threshold and at zero recoil the convexity 
parameter can be seen to go to zero at both ends of the q 2 spectrum because one has 

H ++ — 2H 00 + H —7- at threshold and p ++ — 2p 00 + p — > at zero recoil. 

An interesting exercise is to calculate the LO convexity parameter in the threshold 
region. Neglecting terms of 0(q 2 /m 2 z ), as before, one obtains 



The expression (j65p has the correct threshold behaviour. Keeping q 2 fixed (and small), and 
taking the limit mj — > one has /ij — > 0, A — > 1 and one obtains Cf = —3/2 in agreement 
with Fig. H 

In Figs. [roHT2l we decompose the total differential rate dT/dq 2 in terms of the three 
partial unpolarized transverse (U), longitudinal (L) and scalar (5*) contributions dTjj/dq 2 , 
dTi/dq 2 and dTs/dq 2 , where the three partial rates are defined by the contributions of the 




(65) 
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Figure 10: Differential rate dYjj /dq 2 . Labelling of curves as in Fig. E] 
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Figure 11: Differential rate dTi/dq 2 . Labelling of curves as in Fig. [8] 
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Figure 12: Differential rate dTs/dq 2 . Labelling of curves as in Fig. [8] 

density matrix elements p ++ + p , poo and p tt , respectively. The total rate is then given 

by dT/dq 2 = dTu/dq 2 + dT L /dq 2 + dY s /dq 2 . 

Fig. [10] shows that the transverse rate is weighted toward higher q 2 values, whereas 
the longitudinal rate is more evenly distributed (Fig. [TT|) . The scalar rate is considerably 
smaller and shows a peak close to threshold (Fig. IT2"|) . The peak value is strongly enhanced 
by the radiative corrections. The radiative corrections to the transverse rate are small. 
The radiative corrections to the longitudinal rate can be seen to be quite pronounced close 
to threshold which, in part, is due to the increase of a s due to running. 

In Table [1] we present our numerical results for the integrated total rate and the inte- 
grated partial rates. One can see that the integrated longitudinal rate Tl slightly dominates 
over the integrated transverse rate IV. The scalar rate is quite small and contributes to 
the total rate at the 2.9 % level. The LO total rate is reduced by 5.8 % through mass effects 
where the biggest reduction comes from the longitudinal rate (11.7%). Radiative correc- 
tions increase the LO rates by 6.2 % — 7.8 % except for the scalar rate which is increased 



33 





Born mj 


= 


Born rrii ^ 


0{a s ) 


W*+ — ^ rh 
VV r CU 












r 


2.42 ■ 10" 


-7 


2.28 • 10~ 7 


A A 1fl" 
Z.44 • 1U 


-7 




9.75 • 10" 


-8 


9.35 • 10~ 8 


y.yo • iu 


-8 


T L 


1.45 ■ 10" 


-7 


1.28 ■ 10~ 7 


1.38-10" 


-7 


T s 







6.64 • 10" 9 


6.78- 10" 


-9 


A 

^FB 


n 
u 




U.UZOo 


0.0253 




7* _i. hh 












r 


1.49 ■ 10" 


-5 


1.19 • 10~ 5 


1 Q/1 1fl" 


-5 




6.04 • 10" 


-6 


5.02 • 10~ 6 


O.Oo • 1U 


-6 


T L 


8.87- 10" 


-6 


5.90 ■ 10~ 6 


6.68- 10- 


-6 









1.02 • 10- 6 


1.14-10" 


-6 


■"-FB 


n 

u 




n 

u 


7.39 -10" 


-4 


Z* ->■ cc 












r 


1.16 ■ 10" 


-5 


1.13 • HT 5 


1.20 • 10" 


-5 


IV 


4.69 • 10" 


-6 


4.59 • 10~ 6 


4.84 • 10" 


-6 


r L 


6.89 ■ 10' 


-6 


6.40 ■ 10~ 6 


6.85- 10- 


-6 


r 5 







3.11 ■ 10~ 7 


2.84- 10- 


-7 


^4fb 










5.65-10- 


-4 



Table 1: Integrated rates r,r C /,r i ,r S ' and forward-backward asymmetry A FB for H — > 
W~ + W* + (^r cb), H -r Z + Z*(-^f 66) and H -> Z + Z*(^r 66). All entries are given in 
units of GeV except for A FB . 



only by 2.1 %. We also list the value of the forward-backward asymmetry Afb which, as 
has been discussed before, is a parity-odd effect contributed to by the parity-conserving 
scalar-longitudinal interference term. The forward-backward asymmetry is positive (see 
Eqs. (155]) and ( 16 ip ) and receives its main contribution from the region close to threshold. 
Aps is of the same order of magnitude as Ts/Tl- 

Quark mass effects can be expected to play a larger role in e.g. the decay H — > Z+Z*(-^r 
66). First, the 66 threshold is higher than the cb threshold, and second, the phase space is 
reduced due to the larger mass of the Z boson, i.e. the physical q 2 range becomes smaller. 
An extra bonus is the fact that the decay Z* — > bb is not CKM suppressed. For the 
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differential decay distribution one obtains (sin 2 Ow = 0.23188) 
dT 1 g* 1 \p z \\p\ 1 



dq 2 2 1024vr 3 cos 4 6 W m 2 H ^ (g 2 - m|) 2 + m\Y\ 
\ P++ + poo + p-)(y}H™ L + a)H^ L ) 



x - 
3 



+ 3(l-4) 2 P«(^r + 4^ A )), (66) 



m 



z 

where the gauge boson momentum now is \pz\ = y Mm^j, m 2 z , q 2 ) /2m H , and p w ,m w are 
replaced by pz,raz in the expressions for p mm i in Eq. (|59|) . The electroweak coupling 
coefficients are given by 

g 

Vf = 1 sin 2 6*^, a/ = 1 for u, c, t 

3 

4 

Vf = — 1 + - sin 2 6w, cbf = — 1 for d,s,b. (67) 

In Fig. [13] we provide a plot of the differential q 2 rate for H Z + Z*[— > bb) where we use 
m>z = 91.188 GeV, Tz = 2.495 GeV. Again the differential LO rate shows the appropriate 
threshold behaviour for ^ 0, i.e. the differential rate vanishes at threshold q 2 = 4m 2 . 

The corresponding mj = LO rate shows no apparent vanishing at threshold for the 
same reason as in the corresponding H — > W~W* + case. The differential rate at q 2 = 
and for rnj = is given by 



dT 



= 1 «\ L-^ ^fF^ ^ = 1.120- 10-^GeV-'. (6 f 
g 2=o 2 1024tt 3 cos 4 W 3 m 3 H m 2 z(m 2 z + V%) 2 v 



dg 2 

in agreement with Fig. [131 

As Fig. [13] shows, the NLO rate does not go to zero at threshold. As in the charged 
current case this can be traced to the presence of the NLO chromodynamic Coulomb 
singularity at threshold. One can estimate the finite threshold value of the 0(a s ) rate by 
neglecting terms of 0(q 2 /m 2 z ) in Eq. (I6"6"j) whence one can express the finite threshold value 
in terms of the LO = contribution in Eq. (1681) . One then obtains 



dT_ 

dq 2 



(vj + aj) 16vr 2 dT 
thresh " s (v 2 f + 3a 2 ) 3 ^ dq 2 
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(69) 

g 2 =0 




Figure 13: Differential rate for the three body decay H — > Z + Z*(— > 6 + 6). Labelling of 
curves as in Fig. El 

Note that the contribution proportional to 3a 2 results from the scalar contribution in 
Eq. f )66|) . By a visual inspection of Fig. [131 the approximation can be seen to be quite 
good. Similar to the calculation leading up to Eq. fl65|) one can calculate the LO convexity 
parameter in the threshold region. Neglecting again terms of 0(q 2 /m 2 z ) one finds 

c, = -?i^ (70) 

which is just the limiting case of Eq. ( 165]) for jii = fi 2 '■= Curiously the intricate 
dependence on the electroweak coupling parameters c/ and a/ has dropped out when 
taking the ratio. In the mass-zero case and at q 2 = one has exactly c/ = —3/2. 

In Fig. HHwe show a plot of the q 2 dependence of the convexity parameter Cf. In the 
threshold region the convexity parameter behaves very differently for = and ^ 
(the radiative corrections are quite small). This implies that the polar angle distributions 
are very different for the two cases. In order to illustrate this effect we choose q 2 = 150 GeV 2 
and, in Fig. [151 plot the corresponding cos 9 distribution. At this value of q 2 one is well 
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Figure 14: q 2 dependence of the convexity parameter Cf(q 2 ) for H — > Z + Z*{— > b + b). 
Labelling of curves as in Fig. [8] 

outside of the nonperturbative threshold region. Since the convexity parameter is negative 
(see Eq. (1701). one has a downward open parabola. Mass effects can be seen to be crucial 
for the correct description of the m ^ angular decay distribution which is much flatter 
than the m = distribution. The three curves correspond to convexity parameters of 
c f = -1.154 (LO;m b = 0), c f = -0.388 (LO) and c f = -0.407 (NLO). 

In Table Q] we have listed numerical values for the various integrated partial rates and 
for the asymmetry parameter Aps for both Z* — > bb and for Z* — > cc. Lepton-mass and 
off-shell effects can be seen to be quite large in particular for the bb case. In the bb case, 
mass effects decrease the LO rate by 20.1 % where most of this reduction comes from T^. 
The scalar contribution amounts to 8.6 % of the total contribution. The radiative correc- 
tions increase all four rates by O(10 %). The scalar-longitudinal interference contribution 
sets in only at 0(a s ) since the Born-term contribution to Hot vanishes, i.e. the forward- 
backward asymmetry is proportional to a s and therefore small. This is borne out by the 
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Figure 15: Polar angle distribution for Z*{— > b + b) at q 2 = 150 GeV 2 . Labelling of curves 
as in Fig. M 

tiny numerical value of Afb listed in Table [TJ The numbers in Table [1] for the cc case 
follow a similar pattern, though quark mass and off-shell effects are smaller. 



10 Summary and Conclusions 

We have calculated the NLO QCD corrections to the polarized decay functions in the 
decay of an off-shell and on-shell polarized W + gauge boson into massive quark-antiquark 
pairs VT + (t) — > Q1Q2, keeping the quark masses finite. Using these NLO results for the 
decay process as well as previous results on the NLO corrections to the production process 
t — > b + W + we have studied the NLO corrections to the polar angle decay distribution 
in the cascade decay t — > b + W + followed by W + — > q\q2- We have found that the NLO 
final-state corrections to the decay distribution are somewhat larger than the NLO initial 
state corrections. Altogether we find that the NLO corrections lead to a flatter angular 
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decay distribution W(Q). 

The decay analysis was done in the W + rest frame which has the maximal sensitivity 
to W + polarization effects. Polarization effects of the W + boson will be visible also in 
other reference frames such as the laboratory frame. It is therefore always important to 
retain W + polarization effects in radiative correction calculations (see e.g. Ref. [44J). 

We have presented our results in a general form involving the spin and spin 1 pieces 
of the (VV), (AA), (V A) and (AV) current contributions separately. Our results can thus 
also be applied to on-shell Z decays and off-shell Z* decays (as in Sec. 9) and also to 
extensions of the SM. 

In this paper we have discussed the decays W + — > q± q 2 of positively charged W + bosons. 
The corresponding results for negatively charged bosons W~ — > q~\ q 2 can be obtained from 
the CP invariance of the interaction. One finds [T4] 

H mm (W~ -> q x q 2 \ /ii, /i 2 ; z \\ q 2 ) = H mm (W + -> q x q 2 \ fx 2 , Hi] z \\ qi). (71) 

From the experimental point of view, the leptonic decay of the W boson is the most 
interesting one. In a sequel to this paper we shall calculate the corresponding NLO elec- 
troweak corrections to the decay W /+ (f) — > £ + vt- 

As a further example of much topical interest we have discussed the Higgs decay modes 
H W~ + W* + {^ be) and H ->■ Z + Z\-^ 66, cc) involving the off-shell W* + and 
Z* bosons. We find that lepton-mass and off-shell effects affect the rate and the angular 
decay distributions in these decays in a non-negligible way especially in the vicinity of the 
threshold region. Off-shell effects are also non-negligible for the overall rate. For example, 
the scalar rate T s makes up 8.6 % of the total rate for H — > Z + Z*(—> 66). It would be 
worthwhile to exploit the knowledge about charged and neutral current spectral functions 
in the heavy quark sector which has been accumulating over the last few decades for a 
precision analysis of the rates of the decays H — >■ W~ + W* + (— > be) and H — )■ Z + Z*(— > 
66, cc). 
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A Decay rate terms 



In this Appendix we present analytical expressions for the polarized decay functions intro- 
duced in the main text. For the tree-graph contributions we define logarithmic decay rate 
terms 

ti = In ( l±azttj±vg V ^ = ln ( Wi+iWA \ 

\\ + - n 2 - vAy \i - Hi + fj, 2 - v Ay 

e^Jd^vm, ^ ln ( d + ^-M (A1) 



A*2 y v v/ 1 ! 

and the linear combination £ 3 = t\ + £ 2 . One further has dilogarithmic decay rate terms 
given by 



ij(0) = Li 2 (-, + ) - Li 2 (-,_) + Li 2 J - Li 2 ^^-^ J , (A2) 

SftO) = Li 2 f l-«.-W->A \ + Lb / l-„ 1+fe -v/A \ + / l + ^-^-VA N 
\1 - /ii - /i2 + v A/ \1 - + A*2 + v A/ \1 + Hi - H2 + v Ay 

7T 2 , 1 ^2 / 1 - ^ - fl2 - VX \ + ^ / A \ ^ f 1 - ^ - fl 2 - V\\ 



2 2 \i - /j,i - ^ + V^A/ \2^i/i2y V 1 - A*i - ^2 + V^Ay 
+ 2 ln(2VA*i") ln(2VA*2") - 2 ln(l - ^ + fi 2 + y/X) ln(l + ^ - ^ 2 + \/X), (A3) 



if (0) = Li 2 (//i) - Li 2 (>/7*i2+) - Li^^T^-) - '— 



7T 



2 
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+ l -U 2 ( ^-^l ) + -Li 2 (, 2 ) - 2Li 2 ( VMVE^T \ 

2 2 V(i-^-)V 2 21 2 ^ l-v^r^ ; 

+ In f I— ^ In f 1 Z "~^? r > ) + lnz_ \n{z + - z_), (A4) 
5f(0) = Li 2 (z_) - Li 2 (-z_) - ^ + ln^_ In ' " " ' 



4 \ 1 + z_ J 

- Li 2 + Li2 -i i+ ^f-" n , ( A5) 

I e (0) = Li 2 (^ + ) + Li 2 (V^-) - 2112(7^1) + In 2 z_ = S*(0) (A6) 
where 

z+ = {l + fi! - fi 2 + V\) = zZ 1 . (A7) 



2v^T 

The decay rate terms originating from the loop corrections read 

£ A = 21nA-31n v //H7I 2 -, (A8) 

i B = mg), (A9) 

^ = (^-ln(l-(v^T- v^) 2 ))^ 

+ Li 2 (l - a+) - Li 2 (l - a_) - 2 Re L'^i, /i 2 ) (A10) 

where 

a + = - y= = a_ . (All) 

1 - - n 2 - V A 

The complex function Z/(/xi,// 2 ) is given by 



12 v ' -y/77 ~ / 12 v ~ -vt^ 



where 



+ in ( ln ( VE) (A12) 



m = Li 2 (^) - Li 2 + ?7r .„ (l^f) - r- (A13) 
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and where the velocity parameter v has been defined in Eq. (I33I) . The dilogarithmic and 
double-logarithmic terms in Eq. ( 1A12j) are real whereas L(v) is a complex function with its 



real part explicitly given in Eq. ( 1A13I) . In the limit = ^2 = ^ all dilogarithmic terms and 
the double-logarithmic term in Eq. flA12j) vanish and one remains with the contribution of 



L(v) where v = a/1 — 4/x is the usual velocity of the quarks. Note that the term L'(jii, /i 2 ) 
is a generalization of the equal-mass term (/ii = /i 2 = fi) 

m = Li 2 (^-) - Li 2 (^) + * !„ (1^) - ^ (A14) 

appearing in e + e~ —> tt (see e.g. Ref. [19]). 

B Decay rate terms in the high-energy limit 

In the high-energy or, equivalently, in the mass-zero limit one obtains 

4 -» -ln// 2 , 
£1 — >■ — ln/xi, 

£ 2 ->> -ln/i 2 , 



'3 



— > — In /ii — In // 2 , 



4 -^ln/ii (Bl) 



using the expansion ( )45i) . Further one has 



2;+ -> — — , z- ->■ aT/Ii (B2) 



or, more precisely, <JJ!\Z + — >• 1 — /i 2 . Finally, in the tree-graph case, one obtains 
ij(0) -> -y - ^ln 2 /ii - ^ln/iiln/i 2 - -ln 2 //2, 

7T 2 1 1 

^(0) -> -y - -ln/iiln/ia - -ln 2 /i2, 
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Sf(0) 



S £ (0) ->■ 



/ £ (0) ->■ 



7{(0) 




(B3) 



2 



For the decay rate terms deriving from the loop corrections one has 



-(ln/^i + ln/i 2 ) 



£ B -> In//!— ln/i 2 , 



7r 2 + In 2 /ii + ln/xi ln/i2 + hi 2 /x 2 . 



(B4) 



Finally, one obtains A$ — > 3/4 and A/ — > 3/4. 

C Decay rate terms close to threshold 

Close to threshold where \/\ — >■ one has £ 0) ^2, ^3—^0 while £ 4 — >■ In 4. Note, however, 
that £4 is always multiplied with A and, therefore, does not give any contribution in this 
limit. In order to calculate the dilogarithmic decay rate terms in this limit, one has 
to expand y/X more carefully. To that end we define a small quantity k where k 2 = 
(1 — y/JIi — y/JJ^)- On expanding in k one obtains 




Using the expansion Oil , one can verify that If (0), S*(0), Jf(0), 3{(Q), i*(0) ->■ 0. Finally, 
the decay rate terms that originate from the loop corrections read 

I A -> 2 In A- 31n(V/^i(l - V/*D) > 





(CI) 




t A ->■ 2tt 2 . 



(C2) 
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The term I a appears to be singular at threshold when A — > 0. However, I a is multiplied 
with y/X in Eqs. (HQ]) or £3 in Eq. ( lAlOj) . Therefore, one finds that Aj and As are finite, 



Aj, A s ->• 2vr 2 v ^T(l - (C3) 

Note that the chromodynamic Coulomb singularity at threshold proportional to a s mani- 
fests itself in the overall factor 

N = — S j=N c CfQ 2 ■ (C4) 

7TV A 

D Comparison with spectral function results 

There have been claims and counter-claims in the literature about the correctness of previ- 
ous results on vector and axial- vector spectral functions at 0(a s ). The present calculation 
gives us the opportunity to check on previous results in the literature. According to the 
decomposition 

- gT = -sT + ^f-^f (Di) 

we define the vector and axial-vector spectral functions (H VV ( AA > = H^^^^—g^)) 

H VV(AA) = H VV(AA) _ H VV(AA) (m) 

Following our previous work it is convenient to define the linear combinations (not to be 
confused with the linear combinations H\ and H 2 defined in Sec. 4) 

Hi = \{H™ + H AA ), Hi = \{H™ _ H AA) (D3) 

and, accordingly, for H 1 ' 2 and H^ L . At the Born-term level we obtain 

H\Born) = AN c q 2 (l - m - /x 2 ), H 2 (Born) = 16N c q 2 y/JI^, (D4) 
H l s {Born) = 2N c q 2 (l - ^ - /z 2 - A), H 2 s {Born) = -AN c q 2 ^T 2l (D5) 
H l u+L (Born) = QN c q 2 (l - ^ - fx 2 - A/3), H 2 u+L (Born) = l2N c q 2 ^T 2 . (D6) 
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The NLO corrections read 



N 



4(1 -fi x - fi 2 )A s + 2m(l + nx)l x + 2/i 2 (l + fi 2 )£ 2 



+ ((1 - fix - fi 2 - A)A - 8/xi/x 2 ) 4 - (a*i - fi 2 )\V\£B - 2(1 + /ii + /i 2 - A)VA 



H 2 (a s ) = A^T 2 n\aA s -(3 -fix- 3fi 2 )£ x - (3 - 3fi x - // 2 )£ 2 + 6^" 



(D7) 



and 



H l s (a s ) 



N 



4(1 - fix - fi 2 - A) As - 2fix{fix - fi\ + 16/i 2 - fixfi<i - 4/i 2 Ki 



- 2fi 2 (16fix - Afi\ + fi 2 — fixfi 2 - fij)l 2 — 3 ((1 — fix — fi2 — A)A - 6/ii// 2 ) 
+ 3(/ii - fi 2 )\V\l B + 6(1 - //1 - fi 2 — A + 2//i// 2 )\/a] , 
#f(ai s ) = v^I7^iV[-4A s + (3-/i! -3/X2K1 

+ (3 - 3//i - fi 2 )£ 2 - 6//i/x 2 4 - 3(2 + /ii + /i^v^ . 



(D8) 



Finally, 



N 



H l u+L (a s ) = — [4(3(1- fix- fi 2 ) - \)A s + 2fix(2 + fix + fil-16fi2 + fiil-i2 + 4f4)£x 
+ 2/i 2 (2 - 16/ii + 4/ii + /i 2 + fixfi 2 + fi 2 2 )t 2 - ((1 - fix - fi2 - A)A - 2fixfi 2 ) £3 



+ (fix - ii 2 )XVX£b + 2(1 - 5fii - 5/i 2 - A + Qfiifi 2 )V\ 
Hl +L (a s ) = 3 v ^I7^iV[4A s -(3- /Ul -3/i 2 )£ 1 



(3 - 3^1 - fi 2 )£ 2 - 2fixfi 2 h + (6 - fix ~ A^) 1 /^ 



(D9) 



The normalization factor N has been defined in Eq. f l36l) . The result on Hl J+L (a s ) has 
been listed before in the form 2Hl J+L (a s ) = H u+L (a s ) in Eq. (j4"4"l) . 

When comparing to previous results in the literature we want to remind the reader 
that one uses a different terminology for the spectral function results in the QCD sum rule 
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community. What is called "longitudinal" there is called "scalar" here and what is called 
"transverse" there we call "transverse + longitudinal (U + L)" . 

We find agreement with the results of Ref. [15] which were given in terms of the corre- 
lator functions ImU^ T . These are related to our rate functions by 

H sM = ~— Imn£(s), Hjg(a„) = -— v^iA*2 Imn Z( s )> 

7T 7T 

3]Y 3_/y 

Hh+Ms) = Imn+(s), H 2 u+L {a s ) = y^lmU^s). (D10) 

7T 7T 

We find also agreement with Ref. [16], where the relevant relations are 

lQN c sp y l A {s) = -1^V\(&±H 2 ), 

16N c s p v L /A {s) = JLV\{H$±H* 8 ). (Dll) 



Taking into account the correction mentioned in the note added to Ref. [16] as well as the 
erratum of Ref. [46j, we could not find the obvious mistakes in the integrals J\ and J 2 
mentioned in Ref. [45J. 

We mention that the correlator functions in Ref. jl5l H6] have been obtained by calcu- 
lating the absorptive parts of the pertinent two-loop contributions. The resulting analytical 
expressions for the correlator functions are somewhat simpler than our expressions. The 
mutual agreement was checked numerically. 

E 0(a s ) results in terms of VV, AA, 
VA and AV contributions 

When treating the decay W + — > q\qi we have assumed a SM coupling form for the weak 
decay symbolically written as (V - AY(V - Af = V^V U + N"A V - V>*A U - A»V U . In 
the general case when the relative weight of the vector and axial-vector current is not as 
simple as in the SM charged current transitions (as e.g. in Z — > qq or in SM extensions of 
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the charged current transitions), one wants to be able to avail of the corresponding 0(a s ) 
expressions written in terms of their VV, AA, VA and A V contributions. 

In this Appendix we shall therefore collect all 0(a s ) expressions for the polarized decay 
functions in terms of their VV, AA, VA and A V components. Extending the notation of 
Eq. ([D3]) to 



Hi = \{Hl v + Hi A ), Hi = \{ H yj - H^), 

TT3_i(TjVA TjAV\ IT 4 / TjVA I TjAV\ /T^l \ 

where a is any of U + L, U, L, F, S, tt, tO, Ot, 00, ±± or 1, 2, 3 of Sec. 4, one obtains at LO 

H{(Born) = 2N c q\l - ^ - /i 2 ), H 2 (Born) = AN c q 2 ^T 2 , 

Hl(Born) = 0, H^(Born) = -2N c q 2 V\, 

H^(Born) = 2N c q 2 \, H 2 (Born) = 0, 

Hl t (Born) = 2N c q 2 {l - ^ - fi 2 - A), H 2 t (Born) = -AN c q 2 ^JI^ (E2) 

H^(Born) = -2N c q 2 {^i x - fi 2 ) V\ = H^Born), H 2 (Born) = = H 2 t (Born). 

Using 



ttVV 

H±± 


= Hl + H 2 , 


ttAA 

H±± - 


Hi - Hi 


TjVA 


= ±(H*-iHl) 


ttAV 


= ±(H* + iHl), 


ttVV 


= H\-Hl + {Hl-Hl), 


ttAA ttI 

H oo - H l ~ 




= H tt + H tt> 


Hf t A = 


H tt - H 2 t , 


ttVV 

H to 


= H to + H m 


ttAA 

H to - 


^to ~ H t0 , 


ttVV 


= H ot + H ov 


ttAA 

H ot - 


ttI tt2 

-"ot 11 to J 



(E3) 
one obtains 

Hll{Born) = 2N c q 2 {\ - ^ - fi 2 + 2^T 2 ), 
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H AA (Born) = 2N c q 2 {\ - ^ - ^ 2 - 2^T 2 ), 

Hli(Born) = +2N c q 2 y/\, 

H^(Born) = 2N c q 2 (l - ^ - /x 2 - A + 2 V /^), 

H AA {Born) = 2N c q 2 (l - ^ - /i 2 - A - 2 V /^), 

Hl v {Born) = 2N c q 2 (l - ^ - fx 2 - A - 2 V /^), 

H AA (Born) = 2N c q 2 {\ - ^ - ^2 - A + 2^/JI^), 

Htf) V,AA {Born) = -2N c q 2 {^ - /x 2 )VA = tf 7' AA (i?orn). 



(E4) 



Note that the amplitudes i? 2 ' do not contribute to the parity even pieces of #11^. 
The non-vanishing a s contributions are given by 



H\{a 8 



H 2 (a s ) 



N[2(l fi 2 )A s - 2^(1 + 7/ix - ii 2 )l( 

-v^T(l - 12//i - 5^ - 2/i 2 + 4/i!/i 2 + 
-Ati(6 + 4^! - 7/^)4 + /i 2 (2 + 3^1)^ 
-4/ii/i 2 £ 3 - (1 - 11/ii + // 2 ) VX] , 
^v^II7^[4A 5 + 4/ix/f - 2^(1 + /i! - fi 2 )Si 
-3(1 - ^ - /i 2 )£i - 3(1 - A*i - /i 2 )£ 2 + 3^] , 



# 2 4 (a,) = -7T -4v / AA / + 4(l-3/i 1 -^-2/i 2 + /i2)^ 



N 
~2 

-2(2 - /ii - fij + /i 2 + Ati/^Ko - 8A£ 4 
+4^(1 + 2/i! - /i 2 )£x + 2^(2 + + ii 2 )£ 2 
+ (3 + 14^ - 3//! + 3^) ((1 - v 7 ^) 2 - A* 2 ) 



(E5) 



Hl(a s 



4XA S - 12/^(1 + 7/^ -n 2 )l[ 



-6 v /7ti(l - 12yUi - 5^1 - 2^2 + 4;Ui/i2 + [A)S[ 
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HIM 



-2/ii(20 + 13/ii + ii\ - 24/i 2 + /ii/i 2 + 4/^)4 

+2yU 2 (4 + 12/ii - 4/i 2 - /i 2 - /ilyU2 - ^2)^2 
+ A + /i 2 - (/ii - /i 2 ) 2 ) 4 
-(/ii - ^X^/XCb - 2(3 - 36^i - n\ + 8jj,ifi 2 - 
N^T 2 [l2^l[ - 6^(1 + /xi - n 2 )S{ 
+6//i(l + /i 2 )^i + 6/i 2 (l + //i)£ 2 - 3(3 - //1 — / u 2 )v / A 



(E6) 



H tO/Ot( a s) 



Hto/ot ( 



iV[-2(//i-/i2)VAAj 

+2(/i! - 5/Xi -^1-^2 + ^2 + ^2 + 2/i 2 + /il/i 2 - /i 2 )/^ 

-(3//i - nl — 2u\ - /i 2 - 4/ii/i 2 + 7/i 2 /i 2 + /i 2 + /ii/i 2 Ko - 4(/ii - /i 2 )A£ 4 
+3(/ii - /i 2 ) 2 v / A£i + (4(^i - /i 2 ) + ^2(1 - Mi + ^2) + (Mi - M2) 3 ) v 7 ^ 
- /jlI + /I2 + 2/IX/I2 - hI)\£b =F (j"i - /U 2 )Av / Att - ((1 - v//!!) 2 - /x 2 



x (^5/ii - 8y/JI^fii + 2u\ - 2/x 2 + 2^/7*4^2 - IO/X1/X2 + 2/x 2 



(E7) 

VT^iMa^W + 2(1 - /xi - /x 2 + 3Ati/x 2 )4 - (1 + Mi - Ma)VA4 - A£ B 
±4(/ix - ^v^tt + 3 ((1 - v^I) 2 - M2 ) (1 - 2^/17 - /X! - /x 2 ) ]. (E8) 

The overall normalization factor N has been defined in Eq. ( l36l) . Close to threshold 
\fq I = mi + TTi2 the 0(a s ) results are given by 



H 



H 



+ ^ f-^MiA*2 - ^1^2(16 - 3( v /7i7- y/JI^) (ln/xi - ln/i 2 )) + ©(v 7 !) 



4iV c g 2 {±^y7i77i 2 >0(v / A)} 



3tt 



#0* 



4iV cg 2 ^T-^V^ + C(VA) 



47r 2 a 



Hi, = 4iV c g 2 | ^^7^(VmT-V^)+0(VA)| 
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H 2 0t = -Hi = 4iV c g 2 {-^V^( v ^r- v ^) + 0(v / A)}, 
Hi = -Hi = 4iV c? 2 |v^ 

+ 7^ f^/ii// 2 - V/^i^ (l2 - 3(y7ZI- v /^)(ln/ii - ln// 2 )) + 0(\/A)^ |, 



(E9) 

where, again, identically vanishing contributions are not listed. 
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